Workshop on the CKM Unitarity Triangle, IPPP Durham, April 2003 



m 
o 
o 

(N 
C 

o 



> : 

in . 
o ■ 
m ■ 

\o : 
o . 

m ■ 

o ; 

^ : 

a: 
i 

Oh! 
<D . 
43 ■ 



X 



Charming penguins in 5 — ► 7T7t from QCD light-cone sum rules* 

A Khodjamirian' , Th Mannel, B Melic* 

Institut fur Theoretische Teilchenphysik, Universitat Karlsruhe, D-76128 Karlsruhe, Germany 

We use QCD light-cone sum rules to examine the B — > tttv hadronic matrix element of the current-current operator with 
c quarks in the penguin topology ("charming penguin") as a potential source of the substantial 0(l/mb) effects. Our 
results indicate that charming penguins do not generate sizable nonperturbative effects at finite mj. The same is valid for 
the penguin contractions of the current-current operators with light quarks. The dominant penguin topology effects are 
predicted to be 0(a a ). Still, the nonperturbative effects at finite mb can accumulate to a visible effect that is illustrated 
by calculating the CP-asymmetry in the B® — * ir + iv~ decay. 



1 Introduction 

With the first measurements of the direct CP asym- 
metry in the S-meson system by the BaBar and Belle 
collaborations, charmless two-body hadronic B decays 
have become particularly interesting for constraining 
the CKM matrix elements. In particular, B — > ttt: 
and B — > Kn decays can be used to extract the an- 
gle 7 = arg(V* b ). There are several strategies used 
to determine the 7 angle from such decays, mainly 
based on the isospin and SU(3) relations. Unfortu- 
nately, the theoretical accuracy of these relations is 
limited and it has to be improved by calculating the 
SU(3) breaking effects [[Q. In the direct calculation 
of the relevant hadronic matrix elements one has also 
to resort to approximate methods. Apart from the 
naive factorization, which assumes the vanishing non- 
factorizable interactions, there are methods which try 
to investigate the latter [ [21 [31131- QCD factorization 
approach [ 3 shows that in the rrib — * 00 limit, exclu- 
sive .B-decay amplitudes can be expressed in terms of 
the factorizable part and calculable 0(a s ) nonfactor- 
izable corrections. For phenomenological applications 
it is important to investigate the subleading effects in 
the decay amplitudes suppressed by inverse powers of 
rrifc. Especially interesting are "soft" nonfactorizable 
effects, involving low-virtuality gluons and quarks, not 
necessarily accompanied by an cv s -suppression. Quan- 
titative estimates of nonfactorizable contributions, in- 
cluding the power-suppressed 0(1/ mb) contributions 
can be obtained [|3j using the method of QCD light- 
cone sum rules (LCSR). 

Among the most intriguing effects in charmless B de- 
cays are the so called "charming penguins" . The 
quark pair emitted in the b — > ccd(s) decay propa- 
gates in the environment of the light spectator cloud 



and annihilates to gluons, the latter being absorbed 
in the final charmless state. In this, so called BSS- 
mechanism [ the intermediate cc loop generates an 
imaginary part, contributing to the final-state strong 
rescattering phase. In QCD factorization approach [ 
[3], charming penguins are typically small, being a part 
of the 0(a s ) nonfactorizable correction to the B — > 7T7t 
amplitude. On the other hand, fits of two-body charm- 
less B decays do not exclude substantial 0(l/rrib) non- 
perturbative effects of the charming-penguin type 
Therefore, in [\7\ we have investigated the effects gen- 
erated by c-quark loops in charmless B decays by using 
LCSR. 

2 Charming penguins in LCSR 

The decay amplitude for the B° — > Tr + n~ decay is 
given by the hadronic matrix element (ir + ir~ \H e g \B ) 
of the effective weak Hamiltonian 



H c tf — —= 



c x o\ + c 2 o\ 
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+ ]T Cii^Oi + Cr^ + CsgOsA, (1) 

i=3, ...10 J 

where 0\ = (dT^p) (pT^6) and O p 2 = (pT ^p)(dT>*b) 
are the current-current operators (j> = u,c and V ^ = 
7^(1 — 75)), 03-io are the penguin operators, and CV 7 
and Os g are the electric dipole and chromomagnetic 
dipole operator, respectively. Each operator entering 
Eq. contributes to the B —> 7T7t decay amplitude 
with a number of different contractions of the quark 
lines (topologies). In the discussion we will mainly 
concentrate on the operator 0\. For convenience we 
decompose this operator as 0\ = \0\ + 20%, extract- 
ing the color-octet part &\ = (cr M ^c)(dr^&) . 

The LCSR expression for the B — > tttt hadronic matrix 
element of Of is derived from the procedure described 
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in detail in [ 01 IHJ- ^ ne starts by introducing the 
correlation function: 

= ( P ~k) a F(s l7 s 2 ,P 2 ) + ..., (2) 

where = ujajsd and — imbbjsd are the 

quark currents interpolating pion and B meson, re- 
spectively. Only the color-octet part of the 0\ oper- 
ator contributes at the leading order. The contribu- 
tions of 0\ have to be considered within higher-order 
corrections. 

The LCSR expression for the hadronic matrix element 
A«5S)(B° -> 7T+7T-) = (<K-(p)n+(-q)\d c 2 \B d (p-q)) is 
given by 

f*f B A^XB d ^-K + K-)e- m B' M * = fds 2 e- s ^ 

ml 

xi 1 d Sl e- s i/ M hm S2 Im Sl F( Sl ,s 2 ,P 2 )\ , (3) 

where Mi and M 2 are the Borel parameters in the 
pion and B-meson channels, respectively. The pa- 
rameter Sq (sq) is the effective threshold parameter 
of the perturbative continuum in the pion (B-meson) 
channel. In the sum rule © we take the finite nib 
corrections into account, but neglect numerically very 
small corrections of order s^/m 2 B . The correspond- 
ing sum rule for the hadronic matrix elements of the 
current-current operator with the light quarks 0", 

A { p*\B a d -> 7T+7T-) = (n-Ti+p^Bl) is easily ob- 
tained from LCSR for the 2 operator by putting con- 
sistently m c — > 0. 

We calculate F(si,s 2 ,P 2 ) in © at large spacelike 
si,S2,P 2 employing the operator product expansion 
(OPE) near the light-cone. The corresponding dia- 
grams of 0(a s ) are shown in Fig. ^ They contain 
a c-quark loop, which involves a well known function, 
producing the perturbative imaginary part, due to the 
BSS mechanism The remaining diagrams not 

shown in Fig. with gluons attached to the virtual b 
and d lines, do not contribute to the sum rule because 
their double imaginary parts vanish inside the duality 
regions 0<si < Sq, m 2 <S2< Sq. 

We proceed by investigating the effect of the soft (low- 
virtuality) gluons coupled to the c-quark loop. On- 
shell gluons or light quarks emitted at short distances 
end up in the multiparticle distribution amplitudes 




Figure 1. Diagrams corresponding to the 0(a s ) penguin 
contractions in the correlation function Only the dia- 
grams contributing to the sum rule are shown. The square 
denotes the four-quark operator O 1 ^. 




(a) (b) 

Figure 2. Diagrams containing the quark-antiquark-gluon 
DA's. 

(DA's) of the pion. These contributions are then of 
the higher twist and are suppressed by inverse powers 
of the heavy mass scale with respect to the contri- 
butions of 2-particle quark- antiquark DA's of lower 
twists. Therefore, for our purposes it is sufficient 
to consider diagrams with one "constituent" gluon 
i.e., diagrams involving quark-antiquark-gluon DA's 
of the pion. However, for the given correlation func- 
tion the diagram (Fig. 2a) with one gluon van- 
ishes due to the current conservation in the c-quark 
loop. Nonvanishing terms with the three-particle DA's 
emerge from the diagrams, containing at least one 
hard gluon in addition to the on-shell gluon, Fig. 2b. 
From the studies of the b — > matrix elements of 
01,2, where similar diagrams with an on-shell photon 
and virtual gluon have been calculated [ 9 , it could 
be concluded that the contribution of such diagrams 
is both, a s and 0(l/m 2 ) suppressed with respect to 
the diagrams in Fig. 

Furthermore, from the simple dimension counting 
we find that another possible contribution, shown 
in the diagram Fig. 3a is suppressed by a factor of 
0(l/ml) In(m^). Also, the diagram with two glu- 
ons emitted from the c-quark loop (Fig. 3b) is not 
included in our calculation because it contains DA's 
with multiplicity larger than three and it is suppressed 
with respect to the diagrams in Fig. by at least 
0(l/(m 2 mD). The presence of ln(m^) and m~ 2 in 
the contributions of FigOJi and Fig0t> respectively, 
indicates that at m c — > these terms are divergent 
and by calculating the contribution of 0" from such 
diagrams one will have to consider propagation of the 
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Figure 3. Diagrams corresponding to the multiparticle 
DA's of the pion. and arising from the expansion of the 
c-quark loop in the external ghion field. 





Figure 4. Diagram corresponding to the factorizable 4- 
quark contribution to the correlation function. 



light-quark pair at long distances, i.e. the pion four- 
quark DA diagram [[7|. 

The contributions of four-quark DA's stemming 
from the matrix elements of the type (0 
u{x\)q{x2)q{x3)d(x4)\i{) (xi on the light-cone) are ne- 
glected. But, following we take into account the 
factorizable parts of the 4-quark vacuum-pion ma- 
trix elements, extracting the configurations where one 
quark-antiquark pair forms the quark vacuum conden- 
sate, whereas the other one hadronizes into a twist 2 
and 3 pion DA. Such contributions are enhanced by 
the large parameter /i^ = m\j(m u + m^). In the ap- 
proximation adopted in [0, only two diagrams shown 
in Fig. 01 contribute to the sum rule. Note that the 
quark-condensate diagram in Fig. originates from 
the 4-quark diagram in Fig. OK- Multiparticle con- 
tributions which are factorized in the condensates of 
higher dimension are not taken into account. We also 
neglect the quark-condensate contributions of the type 
(qq)(0 | u(xi)G'j lu (x2)d(x3) \ ir) arising from the dia- 
gram in Fig. Eh after applying the QCD equation of 
motion to the derivatives of G^„. These terms are 
suppressed at least by 0(1/ ml) with respect to the 
diagrams in Fig. 

To summarize, we do not find significant contribu- 
tions involving soft gluons in the OPE of the correla- 
tion function @. The dominant effect which we cal- 
culate arises from the c-quark loop annihilation into 
hard gluons (Fig.Q. In addition, there is the quark- 
condensate contribution (Fig.0J) which we consider as 
a natural upper limit for all neglected contributions of 
multiparticle DA's. 



3 Charming penguins and CP asym- 
metry 

For a numerical estimate of the charming penguin in 
B — > 7T7T decay we calculate the ratio of the sum 
rule © to the factorizable amplitude A^'- ^B^ — > 
7r+Tt-)=im%f w f+ v (0): 



r (Of) ( 53 



7T 7T 



E 

2A^(B° d 



imlUtM {LCSR) 



(4) 



With the parameters taken from [ElCj, and by adding 
linearly the uncertainties caused by the variation of all 
parameters, we get the following for the penguin-loop 
contractions with c and u quarks in the loop: 



= [-(0.29 -r 0.56) - (1.3 1.6)i] • 10~ 2 , 

= [(0.09 ~ 0.21) - (1.6 ± 2A)i] ■ 10~ 2 , (5) 



The charming penguin corrections turn out to be very 
small, not larger than the other nonfactorizable correc- 
tions. However, they appear to produce a noticeable 
effect in the CP asymmetry. Therefore, the CP asym- 
metry appears also to be a good testing ground for 
the influence of the 1 /n%b corrections in the charming 
penguin contributions. 

Penguin contributions influence both the direct and 
the mixing-induced CP violation in B — ► nn. Here 
we concentrate on the direct CP asymmetry in B d — > 
T7 + ir~ decay, which is given as 



dir 
"CP 



lei 2 ) 



(6) 



where £ = e -^(P+i)(\ + R e ^)/(1 + Re-' 1 ' 1 ) and R = 
—P/ (RbT). Here T is the contribution to the B — > tttt 
amplitude proportional to V u bV* d = \V u bV* d \e~ n . It 
contains the tree amplitude, the penguin-loop con- 
tractions of the current-current operators 0" 2 , and 
also V u bV* d proportional penguin operator contrac- 
tions. The remaining contributions, being propor- 
tional to V c bV* d are contained in P. The penguin-loop 
contractions of the current-current operators 0\ 2 rep- 
resent the main contribution to this part. The factor 
Rb = \V ub \\V ud \/(\V cb \\V cd \) is the ratio of the CKM 
matrix elements. 

Both T and P amplitudes have strong phases; there- 
fore we have T = \T\e lST and P = \P\e i5p and the CP 
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Figure 5. Direct CP asymmetry in — > -k + tx~ as a func- 
tion of the CKM angle 7. The upper curve is the result 
obtained for mb — > 00. The dark region is the LCSR re- 
sult, with all uncertainties from the method included (un- 
certainties in the CKM matrix elements are not taken into 
account). The light region shows the deviation from the 
nib — * 00 limit result. 

asymmetry for B® — ► ir + ir can be written as 

dir = -2|E| sm{S P -6 T ) sin 7 

CP l-2|i?| cos{S P -S T ) cos 7 +|i?| 2 ' ( > 

We have to include in T and P all contributions men- 
tioned above. The contributions arising from the pen- 
guin contractions of the operators 0\ and 0" are 
taken as in (J5J • It is easy to extend these results to the 
tree and penguin contributions of the penguin opera- 
tors O3-6. The LCSR result for the gluonic penguin 
contribution of the dipole operator 0$ g is borrowed 
from [ 8 . The electroweak penguin contributions to 
B — > wit are color-suppressed and negligible. 

The hard 0(a s ) corrections to T and P amplitudes are 
known in m b — > 00 limit from QCD factorization [ 3 . 
We have examined the influence of these contributions 
on the phases St and dp. It appeared that they are 
highly suppressed in comparison to the phases emerg- 
ing from the penguin-loop contributions. Therefore, 
we neglect 0(a s ) corrections in 0. 

In Fig. we show function of 7, calculated by 

using the penguin contributions estimated from LCSR 
at finite m,b (dark region) and compare the result to 
the infinite-mass limit that agrees with the QCD fac- 
torization prediction [ 3 j (upper curve). Both results 
are taken at the same scale \j, b ~ m^/l as used in 
LCSR. The prediction shown in the figure is not fi- 
nal, because there are missing annihilation effects, and 
the uncertainty in the CKM matrix elements are also 
not taken into account. However, the figure nicely 
illustrates the size of 0(1/ nib) corrections and the dif- 
ference between the result at the finite and the 
mb — ► 00 result. 



4 Conclusion 

The LCSR estimate presented here for the hadronic 
matrix element of the current-current operator with 
penguin topology involving c and u quarks [0 shows 
that the main contribution to the sum rule stems 
from the 0(a s ) quark loop annihilating to a hard 
gluon, Fig.l. This justifies the generation of the strong 
rescattering phases in B — > 7T7T by the (perturba- 
tive) BSS mechanism. The soft-gluon effects, which 
in the sum rule approach correspond to multiparti- 
cle pion DA's, are suppressed, at least by 0(a s /m\). 
Therefore, we do not find significant nonperturbative 
0(1 /nib) corrections. In mt — > 00 limit our result 
agrees with the QCD factorization prediction for the 
penguin contractions. Since the strong phase is gener- 
ated perturbatively, the CP symmetry in B® — > tt + 7t~ 
is expected to be small. However, at finite mb we show 
that 0(a s /rrib) corrections to a^p accumulate and can 
be noticeable. 
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